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Abstract 

The metal-graphene contact resistance is a technological bottleneck for the realization of 
viable graphene based electronics. We report a useful model to find the gate tunable com¬ 
ponents of this resistance determined by the sequential tunneling of carriers between the 3D- 
metal and 2D-graphene underneath followed by Klein tunneling to the graphene in the channel. 

This model quantifies the intrinsic factors that control that resistance, including the effect of 
unintended chemical doping. Our results agree with experimental results for several metals. 
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field-effect-transistor. 



Figure 1: Sketch of the device considered in this work (a) and the band diagram of a MGOS 
heterostructure (b). Red arrows suggest the current crowding effect near to the contact edge. 


Introduction 

While graphene has emerged as a promising material for future electronic devices thanks to its 
unique electronic properties, the metal-graphene contact resistance ( R c ) remains a limiting factor 
for graphene-based electronic devices. 1 In particular, for high frequency electronics it is an is¬ 
sue very much influencing figures of merit like the maximum frequency of oscillation, the cutoff 
frequency, or the intrinsic gain. 2 Therefore it is neccesary to understand the intrinsic and extrin¬ 
sic factors determining R c , which displays a strong variation depending on the metal contact and 
fabrication procedure details. 3 5 To gain understanding of these factors so that a better control of 
the contact’s technology is feasible, a comprehensive physics based model of R c is an absolute 
requirement. One relevant model was already proposed by Xia et al . 6 to describe the transport in 
metal-graphene junctions as two sequential tunneling process from the metal to graphene over an 
effective transfer length (Lj), followed by injection to the graphene channel (see Fig. |TJi). How¬ 
ever, there is an important ingredient determining R c namely, the transmission from a 3D system 
(metal) to a 2D system (graphene), that has so far not been taken into account properly there in a 
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physics basis. Evidence of the current crowding effect over Lj has been reported by Sundaram et 
al. using photocurrent spectroscopy. 7 

In order to improve the current understanding, we have considered the issue of the carriers 
transmission between materials of different dimensionality. Specifically, we have developed a 
physics-based model where the first process is responsible for the resistance between the metal 
and the graphene underneath (R mg ) and the second process includes the resistance due to a po¬ 
tential step across the junction formed between the graphene under the metal and the graphene 
channel ( R gg ). The total contact resistance is then the series combination of both contributions, 
R c = R mg + R gg , accounting for any current crowding effect near the contact edge. The calculation 
of R mg and R gg are based on the Bardeen Transfer Hamiltonian (BTH) method 8 9 and the Landauer 
approach,® respectively. The BTH method allows us to get information about the matrix elements 
for the transition between 3D-metal and 2D-graphene states and combined with Fermi’s golden 
rule, yields a compact expression for the specific contact resistivity p c . On the other hand, the 
Landauer approach allows to get the conductance of carriers across the potential step between the 
graphene under the metal and the graphene in the channel, where the angular dependence trans¬ 
mission of Dirac fermions and the effective length of the potential have been taken into account. 
To model R c we have considered it as a building block of a FET device, so its value will strongly 
depend on the applied gate voltage. 

Methods 

Electrostatics 

In this paper we start with the graphene electrostatics. We considered a three terminal graphene 
FET (GFET) device controlled by a global back-gate voltage (V g ) as sketched in Fig. |J, al¬ 
though it could be easily adapted to a device with both top- and back-gates as we will show 
later on. We split the electrostatic problem by considering two ID heterostructures, namely, the 
Metal/Graphene/Oxide/Semiconductor (MGOS) and Graphene/Oxide/Semiconductor (GOS) het- 
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erostructures in the contact and channel regions, respectively. In Fig. [1J> the corresponding band 
diagram of the MGOS heterostructure has been shown. In each of these regions we model the gate 
voltage dependence of the graphene Fermi level relative to its Dirac energy, namely AE m and A E g 
for the graphene under the metal and graphene in the channel, respectively. The energy potential 
loops at the encircled interfaces in Fig. [lj) together with the Gauss’s law are considered, resulting 
in Eqs. [TJi-c. Because of the charge transfer between the metal and graphene, a dipole layer of 
size d\ inside the equilibrium separation distance d eq is set up. 11 Also a difference eV between the 
metal and the graphene Fermi level in the contact region, supplied by the drain terminal, has been 
assumed. The work-functions of the metal and graphene are W m and W g , respectively. 


W m = cA c i + W g + A E m — eV, 

(la) 

W g + A E m = eA ox + W sc — e($ s + Vg — V), 

(lb) 

Qsc + Qm + Qg = 0 

(lc) 


In Eq. [lji, the term Aj is the potential drop in the dipole layer which can be expressed as 
Aj = A tr + A c h = —Qm/Cd + A c h, where A tr corresponds to the charge transfer and A c /, to chem¬ 
ical potential interaction describing the short range interaction from the overlap of the metal and 
graphene wavefunctions. 1112 In Eq. |TJ>, the back-gate voltage V g is referred to the source metal 
electrode potential, W sc = X + E g — </>4 is the semiconductor work-function and <j) s is the semicon¬ 
ductor surface potential. In Eq. [lj:, Q m = — QA /r describes the charge per unit area induced in the 
surface metal, Q g ~ 2e/(nfr 2 vj)AE m \ AE m | + Qq is the net charge sheet density within the graphene 
layer 13 plus the charge density due to possible chemical doping 14 ( Qo) and Q sc = C ox A ox describes 
the charge per unit area induced in the semiconductor. Here, Q = £<>/ d\ and C ox = e/T ox describe 
the capacitive coupling to the metal and back gate, respectively. The value of d\ can differ from 
the equilibrium distance d eq (~ 0.3nm) due to the spatial extension of the carbon p z and metal d 
orbitals. The value of A c /, strongly depends on the separation distance d eq and it becomes negli- 
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gible for d eq > 4nm. 11 Combining Eqs. [I] and assuming that <p s saturates at strong inversion and 
acummulation, we get a simple quadratic equation for A E m : 

aAE m \AE m \ T ( C ox T C^) A.E m T cC ox (V g Vd) = 0 (2) 

where a = e 1 j 7r/? 2 v 2 , with vy(~ 1 x 10 8 cm/s) the Fermi velocity, and 

eV D = £?-(W m -W g + eV-eA ch )+(x + eV + e<l> A -W g -f^) (3) 

represents the Dirac gate voltage required to achieve A E m = 0 and defines the back-gate voltage 
value for which p c and the resistance R mg become maximum, as we will see later. 

Because the dipole layer has been modeled as an insulator, the channel region electrostatics 
under the influence of both top- and back-gates can be described in a similar way as presented in 
Eqs. [TJ so the Fermi level shift in the channel (A E g ) can be obtained from: 


aAE g \AE g \ + (C b + C t )AE g + eC b (V b -V gD )=0, (4) 

with 

eVgo = (W m - - eV t ) + (x + eh ~ W g - ^ . (5) 

In the last equations, Cu t \ and Vu t \ are the back (top)- capacitance and gate voltage, respec¬ 
tively. The new Dirac voltage V g o must be understood as the back-gate voltage needed to achieve 
A E g = 0 at a fixed top-gate voltage V t . When there is only a back-gate, like in our experimental 
devices, we can get A E g for the GOS structure simply setting Q —y 0. 

With the electrostatic model given by the above equations, the key quantities A E m and A E g 
could be determined, which in turn are needed to calculate the contact resistance. Fig. [2] shows 
both A E m and A E g at equilibrium (V 0) as a function of the back-gate bias using Palladium 
(Pd) as metal electrode and SiC >2 as oxide with thickness T ox = 90nm. This can be done either 
by solving Eqs. |T| or the simplified Eq. ??, with a very little difference between them. Different 
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kinds of junctions may build-up depending on the back-gate bias, namely pp-type, pn-type, and nn- 
type. Here, we have assumed that Qq is only affecting the graphene channel, and not the graphene 
underneath the metal. The impact of Qo in determining the crossing of A E g with zero, can be 
seen in the figure. To capture the transition between the pp-type and pn-type junction, which was 
observed at V g = V g o ~ 23V, 6 the parameter Qo/e was set to —5.4 x 10 12 cm -2 . Next transistion 
produced between the pn-type and nn-type junction was captured by our model at V g = Vd ~ 46V, 
in accordance with the experiment of Xia et al. The electrical parameters that we have assumed for 
all the simulations presented in this work are shown in Table I. Because of charge transfer between 
the graphene underneath the metal and the graphene in the channel, a potential step of effective 
length A arises at the contact edge. An sketch of that potential step is illustrated in Fig. |2j Once 
we get the electrostatic model, we are now ready to discuss how to model the contact resistance. 



Back-gate bias (V) 


Figure 2: Graphene Fermi level shifts with respect to the Dirac point for different values of Qq 
using Pd as metallic contact. Solid lines: numerical solution of Eq. (1) and symbols: solution of 
Eq. (2). The inset shows the potential step between the graphene underneath the metal and the 
graphene in the channel, with effective length A. 


Resistance R mg and resistivity 

The procedure to model R mg is based on the Transmission Line Method J-^ 17 which in turn requires 
determination of p c , namely: 
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fl mg (A£ m ) = y/^coth(L c /L r )/W, (6) 

where p c (AE m ) = (dJ/dV)~ l |y=o represents the specific contact resistivity, R™ h (250£2/D in 
this work) is the graphene sheet resistance under the metal, Lj = \JPc/R™, is the characteristic 
length over which current injection occurs between the metal and the graphene layer (transfer 
length), and L C {W ) is the length (width) of the contact. Here, p c is calculated by means of the 
BTH method, which allows us to split the metal-graphene system into separate metal and graphene 
subsystems with known Hamiltonians. In the framework of the BTH method, the probability of 
elastic tunneling is calculated using Fermi’s golden rule. This gives a quantitative estimate of the 
coupling between the metal and graphene states, so it is possible to get an analytical formula with 
key parameters for p c as a function of A E m . In the Supplementary data we show how to calculate 
p c from the tunneling current density J using the BTH approach. The resulting compact analytical 
expression for p c as a function of A E m under the metal at V = 0, for a given temperature T can be 
written as: 


Pc(A£' m ) — 


nyh 3 v 2 f exp(2Kd eq ) (£j| - AE m + y 2 E K ) 


\6e 2 


£3(£|l _A£ m ) 2kT ln ( ex P) + 1) - A E n 


(7) 


where y = m/mo, with m and mo the effective electron mass in the metal and dipole layer, 
respectively. The factor K is the electron decay constant in the dipole layer and has the form 
fh 2 +k 2 ,G- 3 where </> ~ W m has been taken as the barrier height and k\\ is the parallel mo¬ 
mentum at the K or K 7 points (i.e., 4k/ 3a). For a typical metal the work-function is W m ~ 5eV 
so K ~ 20nm~ 1 . As a consequence E\\ = h 2 k 2 ^/2m ~ lleV and E K = h 2 K 2 /2m ~ 16eV. From Eq. 
?? we can infer that the maximum value of p c depends exponentially on the equilibrium separa¬ 
tion distance d eq and that the maximum resistivity is located at Vo (A E m = 0). Fig. [3] shows p c 
at T = 300K as a function of the back-gate bias overdrive (V g — Vo) considering different metals. 
After sorting the metals by their peak contact resistivity, it appears that d eq is the main factor con- 
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Figure 3: Specific contact resistivity at room temperature for different metal electrodes centered at 
Vg = Vd- The contact resistivities of Ni and Ti were multiplied by 20 in this plot. 

trolling it, being the Ni contact the best option, followed by Ti. Here we have assumed Si02 as the 
insulator with T ox = 90nm and equal effective masses for every metal. According to Table I, the 
Ni-graphene (Al-graphene) contact has the smallest (largest) equilibrium distance d eq of the metals 
here represented, giving rise to the smallest (largest) value of p c at Vo- The values of p c predicted 
from Eq. ?? are consistent with experimental results reported by Nagashio and Berdebes 18 19 for 
Ni, Ti and Pd. Although the Ni happens to be the best option to get the lowest R c , other effects that 
contribute to the lateral resistance must be considered. As a matter of fact, R c for Ni can become 
comparable to that of Pd, as we will show later. Our model predicts how p c depends on factors like 
the workfunction difference, the equilbrium distance, the chemical interaction potential, the gate 
capacitance and the temperature. 

Resistance R gg 

Next, we model the lateral contact resistance R gg across a potential step with effective length A 
(see inset of Fig. [2]) relying on the Landauer approach. The potential along the transport direction 
x can be described by a simple space-dependent Fermi level shift: 10 
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Table 1: Electrical parameters for selected metal electrodes. They were extracted from previous 
reports, 6 11 except the quantities marked with ” * ” which were considered as fitting parameters to 
match the experimental results. 


Metal 

W m (eV) 

deq ) 

A ch (eV) 

Ni 

5.47 

2.05 

O 

oo 

* 

Ti 

4.65 

2.10 

o 

* 

Pd 

5.67 

3.00 

0.90 

Cu 

5.22 

3.26 

0.99 

pt 

6.13 

3.30 

0.93 

Au 

5.54 

3.31 

0.91 

Ag 

4.92 

3.33 

0.88 

Al 

4.22 

3.41 

0.77 


/ x AE„ — AE n 

AE(x)=AE m + 8 


( 8 ) 


exp (—x/X) + 1 ’ 

where we have considered that the metal electrode cover the left half-plane (x < 0). The type (n 
or p ) and density of carriers in both left and right half-planes are tuned by the back-gate. The im¬ 
portant quantity to be determined is the reflection probability of Dirac fermions across the potential 
step, which has been derived by Cayssol et al. namely: 


R 


step 


sinh^/lK - " 1 " )sinh(^A/c + ) 


(9) 


sinh (n X k + +) sinh (n X k ) 

where the momenta K pa = {AE g — AE m ) /hvp + pk$ + ok ^, with p, <r = ± 1. The longitudi¬ 
nal momentum k x is related to the transversal momentum k y by the phytagorean relationship 


kx ] = sgn (A Ei) yj{AEi/hvff^kj, i = m, g. 

where the positive (negative) sign indicates that the doping is p (n) type. 
By means of the Landauer formula the conductance can be obtained from: 


( 10 ) 


i 2e 2 W f kf 

R- g \AE m ,AE s ) = —-J_ t T s 


Stepdky 


( 11 ) 


2e~W r kF 

'-k F 

where T step = 1 —R ste p is the transmission probability and kp = min (|AE,„|, |AE„|) /hvf. Fig. 
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[4 Y shows the transmission probability of the Dirac fermions across the potential step as a function 
of V g for different incidence angles assuming Pd as the metal. In particular, it indicates the absence 
of backscattering at normal incidence ( k\ = k F or 6 = 0), because of the orthogonality of 
incoming and reflected spinor states. In contrast, the transmission of the bipolar contacts (case pn) 
tends toward zero for incident carriers when d Jt/ 2. 



Figure 4: (a) Transmission probability through a potential step with A = 4nm for different angles 
of incidence. Lateral contact resistance R gg for different values of t\ and h ~ lOOmeV (b) and 
different values of A (c). We have considered here Pd as the metal electrode. 


So far we have not considered the effect of the drain bias (V d ) in defining the contact resistance 
at the drain side ( R c d )• However, for Radio-Frequency (RF) applications, V d is usually placed in 
the saturation region, so its value could be high as compared with V g . In such a case the drain 
and source contact resistances can be very different. Our model for R c is still valid and useful to 
determine R cd in this situation. For this purpose it would be needed to evaluate R c at the effective 
gate voltage V gfiff = V g -V d instead of V g , namely R cd = R c (V gteff ). 


Results and discussion 

Until now, in the description of our model, we have not taken into account any broadening to the 
graphene states in the R c model. To get a more realistic model, an effective broadening describing 
the coupling between the metal and the quasi-bounded graphene states underneath and/or the spa¬ 
tial variations of the graphene-metal distance in the contact surface,^- 1 must be taken into account. 
This effect can be considered upon application of a Gaussian function G\ of width t\ (broadening 


10 















energy). In addition, we have included the random disorder potential in the graphene channel using 
a Gaussian function G 2 of width ^ = hvf \/2niU) where no is the minimum sheet carrier concen¬ 
tration. Then, the two components of R c have to be recalculated as shown in Eqs. S16-17 of the 
Supplementary data. 

In Fig. |4j) we show the effect of t\ on R gg when it varies from 100 to 300 meV with ^ ~ 
lOOmeV (no = 5 x 10 11 cm -2 ). R gg exhibits a main peak corresponding to the minimum DOS in 
the channel (A E g = 0 or equivalently V g — V g o = 0) and another secondary peak corresponding to 
the minimum DOS in the graphene under the metal (A E m = 0 or equivalently V g — V g o ~ 23 V). 
According to the experimental data reported by Xia et al . 6 for Pd as metal electrode, the latter peak 
does not appear in the R c curve, suggesting a large t\ (> 300 meV) value as reflected in Fig. (4jx 

As a complementary information, the dependence of R gg on the effective length A of the poten¬ 
tial step between the metal-doped graphene and the gate-controlled graphene channel is presented 
in Fig. |4j:. For unipolar juntions, R gg is almost independent of A while for the bipolar pn junction 
it moderately increases as A changes from 2 to 128 nm. 

After presentation of the R c model, next is benchmarking it against experimental measurements 
in graphene FETs using the transfer length method (TFM) for metal electrodes such as Palladium 
(Pd), Nickel (Ni) and Titanium (Ti) as shown below. 

In Fig. [5] we have plotted the data reported by Xia et al. considering Pd as metal electrode. 
Here the graphene sheet was transferred to SiC >2 of 90 nm thickness. Our model reveals that R mg 
and R gg play a similar role. The absence of a peak in the experimental R c data at V g ~ 46V suggests 
a large value of t\, as it has previously been discussed. To match the experimental data we have 
assumed Qo/e = —5.4 x 10 12 cm -2 , t\ = 300meV, ^ = lOOmeV and A = lOOnm. Interestingly we 
capture the correct value of the Dirac voltage at V g o ~ 23V and the moderate asymmetry between 
the left and right branches: being R c lower for the left branch because of the much better carrier 
transmission of the unipolar pp junction as compared with the bipolar pn junction (see Fig. |4^). 

Next, we compare with experimental data of GFETs with Ni as metal electrode (Fig. [6]). In this 
case, the back-gated graphene transistors have been fabricated by photolithography on Si wafers 
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Figure 5: Contact resistance and breakdown into its components as a function of the back-gate bias 
overdrive for Pd as metal electrode, where X has been assumed equal to lOOnm. 


covered by 300 nm of thermal oxide. Graphene grown by chemical vapor deposition (supplier 
Bluestone Global Tech) was transferred by the standard PMMA method^ to the substrate and 
patterned using oxygen plasma. Nickel-contacts have been fabricated using sputter deposition and 
lift-off technique. The distance between source and drain contacts was 0.6, 0.9, 1.3, 1.7 and 2.7 
/rm for different devices on the chip to allow extraction of R c by TLM. The channel width was 
10/im. Finally the devices have been encapsulated by 85 nm of AI 2 O 3 deposited by atomic layer 
deposition. After some electrical measurements, we report in Fig. [ 6 jr the comparison between the 
experimental data and the usual model of the source to drain resistance Rt given by : 21 


RCh 

R T = ^L ch +2R c . 


( 12 ) 


Here, the channel sheet resistance R c s ^ has been modeled as R c s % = 


HeJnl + n(V g f 


-1 


, with 


/1 = 1793cm 2 V _1 s _1 and no = 5 x 10 11 cm -2 which were extracted from the experiment, and 
n oc Ais the charge sheet concentration in the graphene channel region. In this case we have 
assumed a possible doping concentration in the graphene channel of Qo/e = —2 x 10 10 cm -2 in 
order to capture the position of the Dirac voltage. Details of the electrostatic behavior of the 
Ni-graphene contact can be found in the Supplementary data. For the quasi-static measurements 
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of resistance shown in Fig. [6^t hysteretic behavior is observed, which is typical for graphene 
FETs. This hysteretic behavior occurs mainly because of charge traps generated by adsorbates, 
typically O 2 /H 2 O redox couples, at the graphene/dielectric interface. 22 23 This effect has not been 
considered in this model. Regarding the contact resistance (Fig. @>). our model gives values 
between 150 and 350 Cl-nm, which are consistent with the experimental values extracted by TLM 
for the gate voltages V g = -20,0 and 20V: R c ~ 220,400 and 220 Q-jitm with correlation coefficient 
R 2 = 0.9894,0.9740 and 0.9754, respectively. The values of X and t\ were determined to be around 
4nm and 300meV, respectively, to get R c values in that range. It is worth metioning that R gg is the 
dominant part of R c , which is in contrast with the Pd contact case analyzed before, where R gg and 
R mg played a similar role. 

A third comparison was carried out for GFETs with Ti as metal electrode with geometrical pa¬ 
rameters: L c h = ljLtm, W = 10/im and T ox = 360nm OS'z'bb). Here graphene synthesized by photo- 
thermal CVD on copper was used to fabricate GFETs. 24 25 Regarding the source-drain resistance, 
the experimental data are shown in Fig. [7ji together with the model prediction on Rj- Similarly to 
the Ni case, we have considered the following electrical parameters: fi = 1805cm 2 /Vs, no = 7 x 
10 n cm -2 as extracted from the experimental data. A chemical doping Qo/e = —4.6 x 10 12 cm -2 
was fed in the model to get the position of the Dirac voltage V g o around 75V in accordance with 
the observation. Details of the electrostatic behavior of the Ti-graphene contact can be found in 
the Supplementary data. Using them together with the parameters given in Table I, our model 
results in the contact resistance shown in Fig. |7]b. The calculated R c agrees well with the extracted 
values from TLM at gate voltages V g = —75,15,45 and 75V: R c ~ 500,500,400 and bOOff-pm 
with correlation coefficient R 2 = 0.9932,0.9915,0.9986 and 0.9997, respectively. The values of 
X and t\ were determined to be around 50nm and 300meV, respectively, to get R c values in the 
mentioned range. Unlike Pd and Ni, R c in Ti-graphene contact exhibits a huge asymmetry between 
left and right branches, being R c lower for the right branch. This asymmetry qualitatively agrees 
with measurements carried out for Ti and reported by Xia et al. 6 
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Figure 6: (a) Comparison of the experimental (dotted lines) and simulated (black solid lines) to¬ 
tal resistance between source and drain for a Ni contacted graphene FET. (b) Predicted R c and 
corresponding breakdown into R mg and R gg components. 



Figure 7: (a) Comparison of the experimental (dotted lines) and simulated (black solid lines) to¬ 
tal resistance between source and drain for a Ti contacted graphene FET. (b) Predicted R c and 
corresponding breakdown into R mg and R gg components. 
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Conclusions 


In summary, we have developed a model of the gate tunable metal-graphene contact resistance. 
First of all we have modeled the behavior of the shift Fermi level A E in both the graphene under¬ 
neath the metal and graphene in the channel. A E becomes zero under the metal at gate voltage 
named as Vo which is controlled by intrinsic electrical parameters like the work function, the ca¬ 
pacitive coupling between the metal and the gate and the value of the chemical interaction. In 
the channel region A E is zero at gate voltage V g o which is strongly determined by the unintended 
chemical doping. Once we get A E in each region, we use a combination of the BTH and the Lan- 
dauer formula to independently determine the contribution of each R c component, namely BTH to 
determine the resistance between metal and the graphene underneath (R mg ) and Landauer formula 
for the resistance between graphene under the metal and the graphene in the channel ( R gg ). Using 
BTH we have found a simple analytical expression for the specific contact resistivity p c which 
elucidates its dependence with the metal-graphene equilibrium distance. Specifically, among the 
metals considered here, Ni and Ti exhibit the smallest value of p c at their respectives Dirac volt¬ 
ages Vo■ However, given the voltage dependence of p c and the different Vo value displayed by each 
metal metal, Cu or Pd could show even a smaller p c than that for Ni or Ti depending on the applied 
gate voltage. The calculation of p c is key to get R mg by means of Transmission Line Method. This 
resistance shows a peak at V g = Vo- On the other hand, the lateral resistance or R gg , in principle, 
exhibits two peaks. One of them at V g o and another at Vo- However when a broadening of the 
graphene states under the metal (t\ in this work) is considered, the latter peak could disappear. 
We have also found that R gg is sensitive to the effective length (A) of the junction potential step, 
specially when a bipolar pn junction builds up. Depending on the metal electrode and the chemical 
doping of the graphene channel the two components of R c could be either similar in magnitude or 
of very different order. In particular for Pd those two components compete, but for Ni and Ti the 
lateral resistance is the dominant component. 

Our model is in agreement with experimental data for several metals under test. In particular, 
we have benchmarked the model against experiments using Pd, Ti, and Ni. The proposed model 
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unveils the interplay between different intrinsic and extrinsic factors in determining the contact 
resistance of graphene-based electronic devices, which should be useful for its optimization. 
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Supporting Information 

Calculation of the specific contact resistivity 

In this section we derive the analytical expression for the specific contact resistivity of the Metal- 
Graphene junction given by Eq. (5) of the main text, relying on the BTH approach. The starting 
point is the expression for the tunneling current 

I = 8S8V £ {r gm f g (E g ) [1 - f m (E m )] - r mg f m (E m )[l - f g (E g )\ } (13) 

g.m 

where both the subscripts g and m label the states in the graphene and metal electrodes with 
energies E g and E m , respectively, gs is the electron spin degeneracy, gy is the valley degeneracy, 
and F gm and T mg refer to the tunneling rates for electrons moving from g —y m and m —> g, respec¬ 
tively. Finally, f g and f m are the Fermi occupation factors for the electrons.The tunneling rates are 
given by the Fermi’s golden rule as 

2 Tt o 

r gm = \Mgm | S {Eg — E m ) = , (14) 

where 
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dS 


(15) 


Mgm — 



_ vr/ 
x m 


d^\ 

dz ) 


are the matrix elements for the transition, with mo the electron mass in the dipole layer. The 
terms v F ? (r,z) and v P /n (r,z) represent the graphene and metal electron wavefunctions, respectively. 
Then, inserting Eq. (??) into Eq. (??), the tunneling current can be expressed as 


Atic 0 

I = Sv-z-Y, \ M m\ L fg(E g ) -fm(E m )\ d(E g -E m ). (16) 

n 8,m 

Considering the graphene with two identical atoms per unit cell, labeled 1 and 2, the wave- 
function for wavevector k can be written in terms of the basis functions ^>jk(j = 1 , 2 ) on each 
atom as T^^z) = X\ (kg)<f>ik s (r,z) +^ 2 (kg) < l > 2 k^(r,z). The basis functions have Bloch form, 
<t>jk g (r, z) = exp (/kg • r) ujk g (r, z ) / Va, where ujk g (r, z ) is a periodic function and A refers to the 
contact area. These periodic functions are localized around the basis atoms (i.e., as 2 p z orbitals) 
of the graphene , and Ujk,, (r,z) is expected to vary only weakly along the radial coordinate r in 
the graphene. Thus, we assume that Ujk g (r,z ) = fjk g (r)g(z) and we approximate the radially- 
dependent term fjk,, (r) as numerical constants f\ and / 2.- 13 The z-dependence has the usual de¬ 
caying form g(z) = V*e~ Kz , where K is the decay constant of the wavefunction in the barrier. The 
decay constant K has the form ^J2m(j)/h 2 + k ^, 13 where (j> ~ W m is the barrier height in the dipolar 
layer and k» is the parallel momentum. For graphene, the latter term is essentially equal to the 
momentum at the K or K’ points (i.e., 4k/ 3a) so that K ~ 20nm _1 for W m ~ 5eV. 

Both Xi (kg) an< l X 2 (kg) h ave well-known values for graphene in a nearest-neighbor tight- 
binding approximation, 26 


Xi 

1 

e Kia/2 

X2 

“ 7 ! 

se ±ia / 2 


where a is the angle of the relative wavevector, the upper sign is for the band extreme at the K 
point of the Brillouin zone and the lower sign is for the K’ point, with s = +1 for the conduction 
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band (CB) and -1 for the valence band (VB). On the other hand, the metal electrons can be modeled 
as free incident and reflected particles for z > d and with a decaying exponential for z<d, namely 


¥ w (r,z) = 


w 


te 


K(z-d) 


e ik m -r 


9 —ik(z—d ) _|_ re ik(z-d) 


z<d 

z>d 


( 18 ) 


where t and r are the amplitudes of the transmitted and reflected waves, respectively. As usual. 


the matching conditions 'Bm (r,z) | z=rf - = v F(r,z)| z=rf + andrng 1 (d^m/dz) \ z=( j- =m~ 1 (d'Vm/dz) \ z= d+ 
have to be fulfilled, resulting in t = 2k-/ (k- + iKm/mo ). Thus, the matrix elements for the transi¬ 
tions of Eq. ?? can be written as 


M. 


gm 


2m 0 k- + iK 


4k 7 K _ . . e 
z 0 (a) 


-Kd 


1 


m 


VUda 


dSe^- km > r , 


(19) 


m 0 


where we have defined 0(a) = X\f\ + Xifi • The integral on the right-hand side of Eq. ?? 
approaches the delta-function <5(k g — k m ) when A —> o°, implying the conservation of in-plane 
momentum k: | M gm | 2 oc| A -1 f dSe l - k ^ kl " )r | 2 —» <5£ km = t>k g) k m - Incorporating Eq. ?? into Eq. 
??, we get the following expression for the current 


8 7te ( h 2 


H \2mo 


4k\/k, 


— Kd 


E 

kg Mm Mz 


0 (a) 


kl+[—K 

\m 0 


[fg (Eg) fm (E m )] 8 (Eg E m ) (>k s ,.k„, ■ 


( 20 ) 


The delta Dirac function guarantees that only energy-conserving tunneling processes are pos¬ 
sible. From the Fig. (| 8 h) we observe that k g = ky + q, with ky constant and thus £ k £ q . The 


function | 0 (a) | 2 is | f\ | 2 Ts/)*/! cos (a), where | f\ | 2 is a constant of order unity assumed to 
have no dependence on k g . 

In deriving Eq. ?? we have incorporated both the graphene and metal dispersion relations, 
namely E g = E g (k g ) = Eo±hvfq and E m = E m (k m ,k z ) = A 2 (/c 2 ; — k z )/2m, which we have sketched 
in Fig. ( 8p) for convenience and k 2 = kl + q 2 — 2qk\, cos a. Considering Eq. ?? in the limit of large 
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Figure 8: a) Diagram of the graphene momentum relative q to the K point, b) Metal and graphene 
dispersion relations. 


A, kg = k m = k, and then the equation for the tunneling current becomes 

Qjt - p / fe2 2 i 

lI.\ & ( a >\ 2o> K>lfs( E s)-U E :.)]S(E s -E m ), (21) 

' V / q,& z 

where we have defined the function co (k z ) = k 2 / (kz + (mtc/mo) 2 ). The discrete sums over q and 
k z are converted to integrals using the recipes £ q —> A / (Ik) 2 Jf daqdq and Y,k_ E/2n f dk z . 
After some algebra, the tunnel current density becomes 


J = 7} (k) /// dadqdk z q | 0(a) \ z co (k z ) [f g (E g ) - f m (E m )] 8 (E g - E m ), (22) 


where 7](k) 
be written as, 


_ 16 e ( ft 2 K^/ke Kd 


h l 2mo 71 


. The energy difference appearing in the delta-function can 


where k* 2 = q 2 + 2q(% — 


£ s -£» = ^« 2 -* z 2 )=°, (23) 

cos a) +kl - kp with % = mvf/h,kp = 2m/h 2 Eo. Using the Dirac 
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delta function properties we can write 8 
becomes 


2m 




2m/h 2 8 ( k z — k*)/ 1 k z + k* and Eq. ?? 


J = rj (k) ^2 ff dadqq | 0(a) | 2 7 —^—-—- ^eV 

I{ ' h 2 JJ y n 2 (k? + (mfc/mo) 2 ) dp 


(24) 


Since we are interested in the specific contact resistivity (i.e V —» 0) we have approximated the 
Fermi levels difference by eV jt, where /i = (Ep g + Ep m )/2. Given that q <C k\\ is fulfilled, we 
can approximate k* 2 ~ kj — kp to find an analytical solution for Eq (??). Thus Eq. (??) can be 
expressed as 


2m 


kl-k}) 


eV^J-. 


J = ^jj dadqq |®(a)| 


Now integrals of the type 


(25) 


df , . 

= I q 


exp[(E — p) /kT] 
kT(l+exp[(E-p)/kT])' 


; dq, 


(26) 


0 0 

where E — p = eV /2 + A E — hv/q, have to be resolved for every cone. Thus, the current density 
takes the form 


J = ^jln(raip[(eV/2 + A£)/ir] + l). (27) 

» (»>’/) 

Finally, an analytical expresssion for the specific contact resistivity p c = (dJ/dV)~ l |y=o as the 
Eq. (5) in the main text is deduced. 

Effect of the broadening on the contact resistance 

The two components of the contact resistance must be modified to take into account the broadening 
of the states in both graphene under the metal and graphene in the channel, namely; 
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R m g(AEm) = / R m l(Ei)Gi(Ei — AE m ;ti)dEi, 


( 28 ) 


R~ g \AE m ,AE g ) = jj RJ g 1 (E h E 2 )G 1 (E l - AE m -,ti)G 2 (E 2 -bE g ,t 2 )dE x dE 2 , (29) 

where R mg and R gg are given by Eqs. (4) and (9) of the main text and the Gauss function is 

G(x — XQ\t) = exp[(x-xo) 2 /t 2 ]/(t^n). 

Ni-Graphene junction 

Fig. [9a] shows the shift of the Fermi level respect the Dirac point for the Ni-graphene junc¬ 
tion. Important values are V g o ~ OV and Vo ~ 125V defining the crossover between unipolar 
pp-j unction/bipolar pn-j unction and bipolar pn-j unction/unipolar nn-j unction, respectively. The 
electrical parameters for this simulation have been mentioned in the main text. On the other hand, 
Fig. [9b] shows the transmission probability of the Dirac fermions across the potential step for 
different incidence angles. 


Ti-Graphene junction 


Fig. |10a| shows the shift of the Fermi level respect the Dirac point for the Ti-graphene junction. 
Important value here is V g o ~ 75V defining the crossover between bipolar pn-j unction/unipolar 
nn-junction. The electrical parameters for this simulation have been mentioned in the main text. 
On the other hand, the Fig. |10b| shows the transmission probability of the Dirac fermions across 
the potential step for different incidence angles. 
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Figure 9: Electrostatic behavior and transmission of the Ni-graphene contact: (a) Graphene Fermi 
level shifts under the contact and in the channel with respect to the Dirac point as a function of 
the back-gate voltage V g . (b) Transmission probability through a potential step with A = 4nm for 
different angles of incidence. Electrical parameters for this metal were taken from Table I of the 
main text. 
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Figure 10: Electrostatic behavior and transmission of the Ti-graphene contact: (a) Graphene Fermi 
level shifts under the contact and in the channel with respect to the Dirac point as a function of 
the back-gate voltage V g . (b) Transmission probability through a potential step with A = 50nm for 
different angles of incidence. Electrical parameters for this metal were taken from Table I of the 
main text. 


22 





























(2) Schwierz, F. Proc. 7FFF2013, 101, 1567-1584. 


(3) Huard, B.; Stander, N.; Sulpizio, J. A.; Goldhaber-Gordon, D. Phys. Rev. B: Condens. Matter 
Mater. Phys. 2008, 78, 121402(R). 

(4) Nagashio, K.; Nishimura, T.; Kita, K.; Toriumi, A. IEEE Int. Electron Devices Meet. 2009, 
5424297. 

(5) Russo, S.; Craciun, M. F.; Yamamoto, M.; Morpurgo, A. F.; Tarucha, S. Physica E 2010, 42, 
677-679. 

(6) Xia, F.; Perebeinos, V.; Lin, Y.-m.; Wu, Y.; Avouris, P. Nat. Nanotechnol. 2011, 6, 179-184. 

(7) Sundaram, R. S.; Steiner, M.; Chiu H-Y.; Engel, M.; Bol, A. A.; Krupke, R.; Burghard, M.; 
Kern, K.; Avouris, P. Nano Lett. 2011, 11, 3833-3837 

(8) Bardeen, J. Phys. Rev. Lett. 1961, 6, 57-59. 

(9) Tersoff, J.; Hamann, D. R. Phys. Rev. B 1985, 31, 805. 

(10) Cayssol, J.; Huard, B.; Goldhaber-Gordon, D. Phys. Rev B: Condens. Matter Mater. Phys. 
2009 79, 075428. 

(11) Khomyakov, P. A.; Giovannetti, G.; Rusu, PC.; Brocks, G.; Van den Brink, J.; Kelly, P. J. 
Phys. Rev. B 2009, 79, 195425. 

(12) Chaves, F. A.; Jimenez, D.; Cummings, A. W.; Roche, S. J. Appl. Phys. 2014, 15, 164513. 

(13) Feenstra, R. M.; Jena, D.; Gu, G. J. Appl. Phys. 2012, 111, 043711 

(14) Ni, Z. H.; Wang, H. M.; Luo, Z. Q.; Wang, Y. Y.; Yu, T.; Wu, Y. H.; Shen, Z. X. J. Raman 
Spectrosc. 2010, 41, 479-483. 

(15) Schroder, D. K. Semiconductor Material and Device Characterization, 3rd Ed.; John Wiley 
and Sons, Inc.: Hoboken, NJ, 2006. 


23 



(16) Leonard, F.; Talin, A. A. Nature Nanotech. 2011 , 6, 773. 


(17) Reeves, G. K.; Harrison, H. B. IEEE Electron Dev. Lett. 1982, 3, 111-113. 

(18) Nagashio, K.; Nishimura, T.; Kita, K.; Toriumi, A. Appl. Phys. Lett. 2010, 97, 143514. 

(19) Berdebes, D.; Low, T.; Sui, Y.; Appenzeller, J.; Lundstrom, M. IEEE Trans. Electron Dev. 
2011, 58, 3925-3932. 

(20) Li, X.; Cai, W.; An, J.; Kim, S.; Nah, J.; Yang, D.; Piner, R.; Velamakanni, A.; Jung, I.; 
Tutuc, E.; Banerjee, S. K.; Colombo, L.; Ruoff, R. S. Science 2009, 324(5932), 1312-1314. 

(21) Kim, S.; Nah, J.; Jo, I.; Shahrjerdi, D.; Colombo, L.; Yao, Z.; Tutuc, E.; Banerjee, S. K. Appl. 
Phys. Lett. 2009, 94, 062107. 

(22) Xu, H.; Chen, Y.; Zhang, J.; Zhang, H. Small 2012, 8(18), 2833-2840. 

(23) Lee, Y. G.; Kang, C. G.; Cho, C.; Kim, Y.; Hwang, H. J.; Lee, B. H. Carbon 2013, 60, 
453-460. 

(24) Riikonen, J.; Kim, W.; Li, C.; Svensk, O.; Arpiainen, S.; Lipsanen, H. Carbon 2013, 62, 
43-50. 

(25) Kim, W.; Riikonen, J.; Li, C.; Lipsanen, H. Nanotechnology 2013, 24, 395202. 

(26) A.H. Castro, F. Guinea, N. M. R. Peres, K. S. Novoselov and K. Geim, Rev. Mod. Phys. 81, 
109-160 (2009) 


24 



